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Using a two-dimensional Hubbard Hamiltonian for the three electronic bands crossing the Fermi 
level in Sr2Ru04 we calculate the band structure and spin susceptibility x(Q5 ^) i^i quantitative agree- 
ment with nuclear magnetic resonance (NMR) and inelastic neutron scattering (INS) experiments. 
The susceptibility has two peaks at = (27r/3, 27r/3) due to the nesting Fermi surface properties 
and at = (O.Gtt, 0) due to the tendency towards ferromagnetism. Applying spin-fluctuation ex- 
change theory as in layered cuprates we determine from x{^i^)i electronic dispersions, and Fermi 
surface topology that superconductivity in Sr2Ru04 consists of triplet pairing. Combining the Fermi 
surface topology and the results for x(Q5^) we can exclude s— and d— wave symmetry for the su- 
perconducting order parameter. Furthermore, within our analysis and approximations we find that 
/-wave symmetry is slightly favored over p-wave symmetry due to the nesting properties of the 
Fermi surface. 
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The novel spin-triplet superconductivity with Tc =1.5K 
observed recently in layered Sr2Ru04 seems to be a new 
example of unconventional superconductivity [||. The 
presence of incommensurate antiferromagnetic and ferro- 
magnetic spin fluctuations confirmed recently by inelas- 
tic neutron scattering (INS) § and NMR ^^O Knight 
shift respectively, suggests a pairing mechanism for 
Cooper-pairs due to spin fluctuations. This is further 
supported by the observed non 5-wave symmetry of the 
order parameter. Likely Sr2Ru04 is another example of 
spin fluctuations induced superconductivity. This makes 
the theoretical investigation of ruthenates very interest- 
ing. NMR and polarized neutron scattering [|j mea- 
surements indicate spin-triplet state Cooper-pairing. In 
analogy to ^He this led theorists to conclude that p-wave 
superconductivity is present However, by fitting the 
specific heat and the ultrasound attenuation Dahm et 
al. doubted p-wave superconductivity and propose an 
/-wave symmetry of the superconducting order parame- 
ter. A similar conclusion was drawn in Ref. jsj. Recently 
it has been reported that also thermal conductivity mea- 
surements are most consistent with /-wave symmetry . 

Clearly, it is important to analyze more definitely the 
origin of superconductivity, triplet pairing and also the 
symmetry of the order parameter on a basis of an elec- 
tronic calculation. This is difficult, since there are three 
Ru^+ i2g bands that cross the Fermi level with ^2/3- 
filling of every band in Sr2Ru04. The hybridization 
between all three bands seems to cause a single Tc. 
All bands cross the Fermi level and hence, the cross- 
susceptibilities are not small and play an important role. 
In view of these facts the previous theoretical analysis 
of the gap symmetries and competition between p and d- 
wave superconductivity [|0-12| must be re-examined and 
it is necessary to determine superconductivity within an 
electronic theory and to derive the symmetry of the order 



parameter from general arguments. 

In this letter we present an electronic theory which 
takes into account the hybridization between all three 
bands. We calculate the Fermi surface (FS), energy dis- 
persion and the spin susceptibility x including all cross- 
susceptibilities. Then, we analyze the pairing interaction 
mediated by the spin fluctuations exchange in Sr2Ru04. 
Analyzing experimental results for the ^^O Knight shift 
and INS data as well as the FS observed by Angle- 
Resolved-Photoemission-Spectroscopy (ARPES) [|l^ we 
obtain values for the hopping integrals and effective 
Coulomb repulsion U . Taking this as an input into the 
pairing interaction we analyze the p-, d- and /-wave 
superconducting gap symmetries. The delicate com- 
petition between weak ferromagnetic spin fluctuations 
and relatively strong incommensurate antiferromagnetic 
spin fluctuations due to nesting of the FS cause triplet 
Cooper-pairing. We get singulet da;2_^2-wave symmetry 
is energetically less favorable. 

We start from the two-dimensional three-band Hub- 
bard Hamiltonian 
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where ai^^aa is the Fourier transform of the annihilation 
operator for the da orbital electrons {a = xy^ yz^ zx) and 
Ua is an effective on-site Coulomb repulsion. t\^a denotes 
the energy dispersions of the tight-bindings bands calcu- 
lated as follows: t\^a = — — '^tx cos kx — 2ty cos ky + 
At' co^kxco^ky. In accordance with experimental mea- 
surements of the Fermi surface and energy dispersions 
we choose the values for the parameter set {eo^tx^ty^f) 
as (0.5, 0.42, 0.44, 0.14), (0.23, 0.31, 0.055, 0.01), and 
(0.24, 0.045, 0.31, 0.01)eV for dxy-, d^x-, and d^^-orbitals 
[13|. The analysis of de Haas- van Alphen experiments 
[14 1 shows a substantial hybridization between xz- and 
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FIG. 1. Calculated susceptibility Re Xo^(Q5^ — 0) ob- 
tained from electronic calculations using the hybridized 
bands. The wave vectors = ^) and ^ (0.27r, 0) 
reflect nesting transferred from a and f3 bands and the orig- 
inal tendency of the 7 band towards ferromagnetism, respec- 
tively. These wave vectors play the most important role for 
Cooper-pairing. The inset shows results for the energies of 
the hybridized bands (i — a, (3^ 7). 
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FIG. 2. Results are shown for the static susceptibility at 
T = 50K. The solid curve refers to the total susceptibility 
within RPA resulting from the partial susceptibilities Re Xo 
{i — a\ P' ^ '^') shown by the dashed curves. The Xo refer to 
the diagonal elements of the diagonalized matrix [xq]- Note, 
the pairing wave vectors (Q^, q^) are nearly the same as in 
Fig. 1. Note, the smallness of Xo cind Xo compared to 



yz- orbit als about t_\_ = 0.1 eV, but not with the xy- 
orbital [|l^. However, the observation of a single Tc im- 
plies the coupling between all three bands. Therefore, we 
choose a weak hybridization thyb = O.OleV (hybridization 
between xy- and xz-^yz-oihitals) <C t±. Note, even such 
a weak hybridization transfers the nesting properties to 
the xy-orbital. In the inset of Fig. 1 we show the resul- 
tant energy dispersions of the obtained hole-like a-band 
and electron-like P- and 7-bands after hybridization. Due 
to the small value of hybridization between xy and yz^ 
xz orbitals the dispersion curves and resulting FS (see 
also Fig. 4) look quite similar to the non-hybridized ones 
[|l6|. However, the importance of hybridization between 
these orbitals for the spin susceptibility, x{^^^)^ will be 
immediately seen from the analysis of the latter. The 
susceptibility is given by: 
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where /(e) is the Fermi function and is the energy 
dispersion of the a, /3, and 7 band 0. 

In Fig. 1 we show the momentum dependence of the 
real part of Xo^- While hybridization between bands does 
not affect much the energy dispersion, it changes signifi- 
cantly the susceptibility of the 7-band. In particular, the 
nesting properties of xz- and ?/z-orbitals reflected by the 
peak at = (27r/3, 27r/3) in Xo^ are caused by the hy- 
bridization between xz^ yz and xy bands. Note, without 
taking into account the hybridization one would not get 
the peak at in the 7-band, but only a broad hump as 
discussed earlier [^. The small peak at ^ (0.27r, 0) is 
due to the original tendency towards ferromagnetism of 
the x?/-band and is not affected by the hybridization. 

Our results for x{^.^^) will have important conse- 
quences and were obtained in contrast to previous studies 
from observed band dispersions. The susceptibility ma- 
trix [x*-^] is calculated where i, j refers to the hybridized 



bands. As already mentioned the non-diagonalized sus- 
ceptibilities are not small and thus cannot be neglected. 
After diagonalizing the matrix [Xq] we get within RPA 
with an effective U (q) for x* 
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where now Xo(q,^) = Ei/Xo(q.^)- Here, Xo(q.^) 
{i^ = a' ^ P\ and 7') are the diagonal elements of the 
diagonalized matrix [xo ]• This obtained susceptibility 
characterizes the normal state magnetic properties of 
Sr2Ru04. Its spin fluctuations are given by x{^.^^) with 
peaks at Qi and q^. These wave vectors are important for 
determining the symmetry of the superconducting order 
parameter. 

In Fig. 2 we present the results for x{^^^) obtained 
from Eq. (^ and the susceptibilities Xo (q^ ^) obtained 
after diagonalization of [xq (q, u;)] . Here, we approximate 
U (q) hy U = 0.345eV which gives agreement with INS 
[|l8| . Remarkably, the peak in Re x at remains nearly 
the same as for Xo^- The peak at q^ is also present, but 
slightly shifted to a larger value. Clearly, Xo much 
larger than Xo' and Xo • Our calculations also have 
shown that cross-susceptibilities Xo 7^ j) cannot be 
neglected. 

In Fig. ||(a) we compare our calculation of the tem- 
perature dependence of the uniform spin susceptibility 
x(0,0) which is measured by the ^^O Knight shift 
and in Fig. ^(b) we compare Im x(Qi,^) with INS data 

. For the calculation of x(0, 0) we approximate U (q) by 
U{0) = 0.177eV [|l9| which gives agreement with Knight 
shift measurements and is also taken in previous calcula- 
tions. These comparisons shed light on the validity of our 
results for x(q, Note, we also take into account that 
there are four electrons per three t2^-bands that would 
give every Xo an additional weight 4/3. Our results are in 
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FIG. 3. (a) Calculated temperature dependence of the 
uniform spin susceptibility with Uo=0.177eV is compared 
with the ^^O Knight shift measurements. The peak is due 
to thermal activation involving 7 and a, [3 bands, (b) Calcu- 
lated frequency dependence of Im x(Q^,a;) compared to INS 
data using Uq. = 0.345eV. 



fair agreement with experiment that shows a tendency to- 
wards ferromagnetism [Q. The maximum in Xi?PA(0, 0) 
at about 25K results from thermally activated changes in 
the populations of the bands near Ep- In Fig. 3(b) we 
compare our results with INS data. In this case we must 
take Uq. = 0.345eV in order to fit xIq? ^) to the peak 
position and height at 00= 6meV as observed in INS. 

While an uncertainty in the INS data (shown in Fig. 
^(b)) is present, our results for x(Q7^) should be a use- 
ful basis for further calculations. The antiferromagnetic 
spin excitations result in incommensurate antiferromag- 
netic Ru-spin alignment at distances larger than near- 
est neighbors. Hence, if Cooper-pairing involves near- 
est neighboring Ru spins also incommensurate antiferro- 
magnetic fluctuations will cause triplet-pairing because 
neighboring Ru spins see also partly a ferromagnetic en- 
vironment. Note, x(Q7^) controls the symmetry of the 
superconducting order parameter. 

For the analysis of superconductivity in Sr2Ru04 we 
take into account that experiment observes non-5-wave 
symmetry of the order parameter which strongly suggests 
spin-fluctuation-mediated Cooper pairing. Assuming the 
spin-fluctuation-induced pairing it is possible to analyze 
the symmetry of the superconducting state from the gap 
equation and our calculated results for x{^.->^) with the 
pronounced wave vectors at and q^. We get for the 
gap equation: 
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where El^ 



\fe\^ ^ + are the energy dispersions of 



the bands and the pairing potential F^-^-^(k, k') is 
different for singulet (a = 0) and triplet (a = 1) Cooper 
pairing. The eigenvalue analysis of Eq. (4) will yield 
the symmetry with lowest energy. The 7-band plays the 
most important role. 

For the determination of the pairing symmetry we fol- 
low the analysis by Anderson and Brinkmann for ^He 
and by Scalapino for the cuprates [|2^ and use the calcu- 
lated FS and spin susceptibility for Sr2Ru04. For triplet 
pairing the effective pairing interaction is (/7 = /7(Qi)) 
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and for singulet pairing 
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respectively It is important to note that for Ei? 

the gap function can be expanded into spherical har- 
monics corresponding to the angular momentum / = 1, 
2, 3, ... and no mixture of A^ belonging to different 
symmetry representations can be present if a single Tc 
is observed. Therefore, we can exclude the (p + (i)-wave 
superconducting state. Using appropriate symmetry rep- 
resentations 1^ we discuss the solutions of Eq. (^) for the 
p, (i, and /-wave symmetries of the order parameter: 

Ap(k) = Aoz(sin kx -\- i sin ky)^ (7) 
Arf(k) = Ao(cos kx - cos ky), (8) 
A/ (k) = Aqz (cos kx — cos ky ) (sin kx -\- i sin ky) . (9) 

Note, the largest eigenvalue in Eq. (||) will yield the su- 
perconducting symmetry of A^ in Sr2Ru04. Solving Eq. 
(4) in the first BZ down to 5K, we find /-wave symme- 
try slightly favored. As expected p- and /-wave sym- 
metry Cooper-pairing are close in energy (A/ = 0.76 > 
Ap = 0.51). A more complete analysis taking into ac- 
count also the coupling between Ru02 planes and inter- 
band U might yield a definite answer [|2j. Note, to ob- 
tain a combined energy gain from the antiferromagnetism 
and Cooper-pairing one expects an order parameter with 
nodes in the RUO2 plane and possibly also with respect 
to the c-direction. 

The solutions of Eq. (^) can be characterized by Fig. 4 
where we present our results for the Fermi surface, wave 
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FIG. 4. Symmetry analysis of the order parameter for the 
triplet pairing in the first BZ. a, /3, and 7 denote the FS of the 
corresponding hybridized bands. The wave vectors and 
are the pronounced wave vectors resulting from the suscep- 
tibility shown in Fig. 2. These determine the symmetry of 
the order parameter. Also for /-wave symmetry the nodes of 
the real part of the order parameter are shown (dashed lines) 
and the regions + ( - ) where the /-wave superconducting gap 
is positive (negative). Note, for the real part of the p-wave 
order parameter the node occurs along kx=0. 



vectors and and symmetry of the order parame- 
ter in Sr2Ru04. The areas with A/ > and A/ < 
are denoted by (+) and (-), respectively. In a good ap- 



proximation we linearize Eq. 



in A/, i.e. 



and safely put tanh(e^//2/cBT) = 1. Therefore, the main 
contribution to the pairing comes from the Fermi level. 
Note, determining A^ for the Ru02-planes it is sufficient 
to take into account only the 7-band, since only this band 
has a dispersion in the plane. The minus sign in Eq. (^) is 
cancelled for triplet pairing (see Eq. (|^)). Furthermore, 
the summation over in the first BZ is dominated by 
the contributions due to and the one due to the back- 
ground and q^. Thus, we obtain approximately for the 
7-band contribution (/ = / or p) 



(10) 



where the sum is over all contributions due to and 
q^. As can be seen from Fig. 4 in the case of /-wave 
symmetry the wave vector q^ in Eq. (|l^) bridges the 
same number of portions of the FS with opposite and 
equal sign. Therefore, the second term in Eq. ( p!Q| ) is 
approximately zero for triplet pairing. We see from Fig. 
4 that Qi bridges portions of the FS with equal signs of 
the superconducting order parameter. Thus, a solution 
of Eq. (10) for A/ is indeed possible. 

In the case of p-wave pairing the real part of the order 
parameter has a node only along =0 in the k^^ky- 
plane. Then, using the corresponding Fig. 4 the wave 
vectors bridge portion of the FS where ReA^ has the 



same or opposite sign. Regarding the q^ contributions 
the situation is similar as in the case of the /-wave sym- 
metry. Hence, we expect for the eigenvalues < A/ as 
in the result of the algebraic solution of Eq. (4). Note, 
for increasing nesting Fig. 4 also suggests that /-wave 
symmetry is favored more than p-wave. The eigenvalue 
analysis of the possible solutions A/ and Ap + iA/ should 
increasingly rule out the latter for stronger nesting 

Also using similar arguments we can rule out singulet 
pairing on the basis of Eq. (^). In particular, assuming 
(i^2_^2 -symmetry for Sr2Ru04 we get a change of sign 
of the order parameter upon crossing the diagonals of 
the BZ. According to Eq. (Q) wave vectors around 
connecting areas (+) and (-) contribute constructively to 
the pairing. Contributions due to q^ and the background 
connecting the same sign areas subtract from the pairing 
(see Fig. 4 with nodes at the diagonals for illustration). 
Therefore, we get that the four contributions due to q^ 
and the background cancel the pair-building contribution 
due to Qi. As a consequence we obtain no (ia,2_^2-wave 
symmetry. Note, this is in contrast to the cuprates where 
the cancelling contributions due to q^ and the small back- 
ground are negligible. For the da^^y^y^^^^^y where the 
nodes are along (7r,0) and (0,7r) directions we can argue 
similarly and thus exclude this symmetry. 

Thus, as a result of the topology of the FS and the 
spin susceptibility we get for p- and /-wave the strongest 
pairing and can definitely exclude d-wave pairing. In our 
approximation we find that /-wave symmetry pairing is 
slightly favored over p-wave symmetry in Sr2Ru04. 

In view of our Fig. 4 we also remark that while Re 
Af exhibits three line nodes, that can be seen by phase 
sensitive experiments, |Ajp shows nodes only along the 
diagonals, as recently found in measurements of ultra- 
sound attenuation below Tc ||2^. However, note in view 
of the low eigenvalues for p and /-wave symmetries and 
approximations used we cannot definitely conclude that 
/-wave is favored over p-wave. 

In summary, we show that hybridization between all 
three bands is important and transfers the nesting prop- 
erties of xz- and yz-oihitals to the 7 band in Sr2Ru04. 
Taking into account all cross-susceptibilities we success- 
fully explain the ^^O Knight shift and INS data. Most 
importantly, we calculate x{^^^) show on the ba- 
sis of the Fermi surface topology and the calculated spin 
susceptibility x{^.^^) that triplet pairing is present in 
Sr2Ru04. Our analysis seems quite general and also 
predicts d-waxe symmetry for electron- and hole-doped 
cuprates j27|. To decide whether p- or /-wave symme- 
try pairing is present one needs to perform more complete 
calculations including coupling between RUO2 planes, for 
example. If the interplane coupling involves also nesting, 
then corresponding nodes are expected. 
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